Size-dependent melting of spherical copper nanoparticles 
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Abstract 

We report size-dependent melting of spherical copper nanoparticles embedded into silica matrix. 
Based on the temperature dependence of the surface plasmon resonance energy and its width we 
observe two distinct melting regimes. For particles smaller than 20 nm the absorption spectrum 
changes monotonically with the temperature, and this allows us to assume the gradual solid-liquid 
phase transition (melting) of the nanoparticles or existence of superheated solid nanoparticles. In 
contrast, for nanoparticles larger than 20 nm, we observe a jump-like increase of the bandwidth and 
non-monotonic dependence of surface plasmon energy at the temperatures below the bulk melting 
point. This indicates that the melting of large nanoparticles is a first-order phase transition similar 
to the melting of bulk copper. 



PACS numbers: 78.66. Sq, 78.68. +m, 64.70.Dv 



I. INTRODUCTION 



Optical properties of copper nanoparticles are quite remarkable because the energy of the 
dipolar mode of surface collective electron plasma oscillations (surface plasmon resonance 
or SPR) coincides with the onset of interband transition. Therefore, optical spectroscopy 
gives an opportunity to study the particle-size dependence of both valence and conduction 
electrons. The intrinsic size effect in metal nanoparticles, caused by size and interface 
damping of the SPR 1,2 , is revealed experimentally by two prominent effects: a red shift of 
the surface plasmon band and its broadening. Scattering of the conduction electrons by the 
surface boundaries of a nanoparticle causes a dephasing of collective excitations, i.e. SPR, 
and is well described by the expression 1 

r = r 00 + A^, (i) 

r 

where T is width of plasmon absorption band, = l/r^ is the size- independent damping 
constant which is related to the lifetimes of electron-electron, electron-phonon, and electron- 
defect scattering, Vqo is the Fermi velocity in bulk material ( Voo = 1-57 x 10 8 cm/sec in 
bulk copper), r is the radius of nanoparticle, and A is a theory-dependent parameter that 
includes details of the scattering process (e.g. isotropic or diffuse scattering^ 2 - 1 ^) . Note, 
that Eq.(TjQ) describes the width of SPR band only when the onset of interband transitions 
has a higher energy than that of SPR energy, e.g. for silver nanoparticles. The confinement 
of conduction electrons in a small volume of a nanoparticle causes a blue shift of SPR band 
when the size of the nanoparticle is decreased, accompanied by a spill out effect beyond 
the physical boundary of the nanoparticle 1 . The particle-size dependence of the interband 
transition edge is caused mostly by the surface stress-induced changes in the lattice constant 
and corresponding changes in the electronic band structure^. This leads to an additional 
blue shift of SPR band with decreasing particle size. 

Melting of crystalline materials is characterized by lost of the crystal order and following 
by a transformation to a liquid phase. As it has recently been shown, melting of small 
nanometer-size particles is quite different process than that of bulk materials (see, e.g. Refs. 
QJfjQJsl)- As the size of a particle decreases, the surface to volume ratio is increased inversely 
proportional to the radius of the particle. As a result, the surface atoms become more prone 
to detach from their positions and diffuse on the surface. Surface diffusion is followed by a 
partial melting of the nanoparticle mainly within a thin layer close to the surface, which is 
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called surface melting or premelting 9 -. The surface melting generally occurs at temperatures 
much lower than the bulk melting but the details depend on the parameters such as size^, 
morphology^, defects^, strain^, and crystallographic orientation of the surface^. Finally, 
similarly to the surface melting, the interior or volume melting of a nanoparticle takes place 
at melting temperature (T m ) lower than the bulk melting point. The decrease of T m with 
decreasing particle size is a well known phenomenon that has been described in terms of 
phenomenological thermodynamic models^, and, more recently, simulated using molecular 
dynamics*^ 1 ^ 1 ^. Experimental evidence of the decrease of T m has been demonstrated by a 
variety of techniques^ 1 ^ 1 ^ 1 ^. 

We have shown earlier— that absorption spectra of Cu/Si02 nanocomposite reveal a clear 
SPR peak located around 560.8 nm (2.210 eV) - 581.1 nm (2.133 eV) depending on the size 
of Cu nanoparticles. The SPR peak position is in good agreement with that of predicted 
by the Mie theory. As we reported in the previous paper 19 , the dependence of the width 
of plasmon peak on the nanoparticle diameter is perfectly described by the well known 
1/r law (expression (OQ)). The fitting of our experimental dependence T(d) obtained at the 
temperature of 293 K gives the following values of parameters (CQ): Poo,293 = 0.087 eV and 
A 293 = 0.107. 

Here we report the size dependent melting transition in spherical copper nanoparticles in 
a broad size range 5-65 nm by means of the absorption spectroscopy, which is, to our best 
knowledge, has not been reported so far. We noted already that the melting temperature 
of metals decreases when the size of the nanoparticles decreases. Since metal nanoparticles 
have been widely used in electronic and optoelectronic industry, the determination of the 
maximum allowed temperature for their exploitation is apparently a very important tech- 
nological problem. This information is obtained from experiments of metal nanoparticles 
melting. Moreover, the optical spectroscopy is a nondestructive tool for determination of the 
extremal characteristics of exploitation regimes of nanometal composites based electronic de- 
vices. Besides the technological interest, the melting temperature is of fundamental scientific 
interest. In most reports the melting of nanoparticles has been detected by high resolution 
electron microscopy (HRTEM) and electron diffraction techniques. However, both tech- 
niques use the samples which are sufficiently heated by the electron beam. Therefore, there 
is no possibility to control the real temperature of the sample. On the other hand, the ab- 
sorption spectroscopy gives such a possibility, i.e. the temperature determined in our work is 
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the real temperatures of the sample. The study of the temperature dependences of the SPR 
energy and bandwidth in copper nanoparticles in broad size range have allowed us to reveal 
the two principally different regimes of melting of small and large copper nanoparticles. 



II. FABRICATION OF COPPER NANOPARTICLES. EXPERIMENTAL PROCE- 
DURES 

Cu nanoparticles in a SiC>2 matrix were synthesized using a procedure described in Ref. 



20l . Porous silica matrices (obtained by a transformation to monolithic glasses) were pro- 
duced by the conventional sol-gel technique based on hydrolysis of tetraethoxysilane (TEOS). 
We modified the procedure by introducing a doping followed by a chemical transformation 
of the dopants in air or controlled gaseous medium. A precursor sol was prepared by mixing 
of TEOS, water and alcohol, with the acid catalysts HN03 or HC1. Silica powder with 
particle size about 5-15 nm (the specific surface area is 380 ± 30 m 2 /g) was added into 
the sol followed by ultrasonication in order to prevent a large volume contraction during 
drying. The next gelation step resulted in formation of gels of desired shape. Gels were 
dried at room temperature under humidity control, and porous materials (xerogels) were 
obtained. Annealing in air at the temperature of 600 °C during 1 h allowed us to control the 
porosity of SiC>2 matrices. The heating process resulted in visible variations of the density 
and specific surface area. Copper doping was performed by impregnation of xerogels in 
Cu(NO"3)2 alcohol solution during 24 h. Then, the impregnated samples have been dried in 
air at 40 °C during 24 h. Further processing of the Cu-doped xerogels was done in three 
different ways: (1) Annealing in air with gradual increase of the temperature from 20 °C 
to 1200 °C (annealing time at 1200 °C: 5 min). (2) Initial annealing in air with gradual 
increase of the temperature from 20 °C to 1200 °C (annealing time at 1200 °C: 5 min), then 
annealing in the atmosphere of molecular hydrogen at the temperature of 800 °C during 1 
h. (3) Annealing in H2 with gradual increase of the temperature from 20 °C to 1200 °C 
(annealing time at 1200 °C: 5 min). As we have shown in the previous work^, the annealing 
in air leads to formation of both the copper nanoparticles and the copper oxide (CU2O) 
nanoparticles. An annealing in hydrogen results in reduction of Cu(I,II) to Cu(0) that is 
aggregated in the form of Cu nanoparticles^. Glass samples were polished up to thickness 
about 1 mm for optical measurements. The samples obtained by annealing in air have low 
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optical density, and are non-colored or slightly yellow (labeled hereafter as samples A); ones 
obtained at successive annealing in air and hydrogen have higher optical density, and are 
light-pink, orange or red colored (samples AH); ones obtained at annealing in hydrogen have 
the highest optical density, and have red color (samples H). 

A tungsten-halogen incandescent lamp was used as a light source for the absorption 
measurements. The single spectrometer MDR-3 has been used for in situ measurements of 
absorption spectra. Absorption spectra of Cu/Si02 nanocomposites containing Cu nanopar- 
ticles in size range 5-65 nm have been measured in the temperature region of 77 - 910 K. The 
measurements at 77 K were performed for the samples dipped in the liquid nitrogen, and 
measurements at the room and higher temperatures were performed in air. The samples were 
placed into an open furnace during the absorption measurements at high temperatures. The 
thickness of the samples was controlled before and after heating and remained unchanged. 
The measured absorption spectra were unchanged after several cycles of heating and cooling, 
therefore, we conclude, that during the heating of samples performed at absorption mea- 
surements the copper particle diffusion in the matrix is absent. Therefore, the filling factor 
of nanoparticles in matrix and inter-particle distances did not change. So, particle diffusion 
can not be the cause of the effects observed in the temperature dependences of SPR band- 
width and energy considered in section HVl Transmission electron microscopy of studied 
Cu/Si02 nanocomposites was performed by JEOL JEM-2000EX electron microscope at a 
room temperature. 

III. TRANSMISSION ELECTRON MICROSCOPY OF Cu/Si0 2 NANOCOMPOS- 
ITES 

Some TEM images of the nanoparticles synthesized at the different technology conditions 
are shown in fig. [TJ The images prove the formation of spherical Cu nanoparticles of different 
sizes in the matrix depending on the conditions of annealing. The copper nanoparticles of 
small size (diameter d is in the range of 2 - 15 nm) were grown in the samples annealed in 
air (samples A). The medium-sized nanoparticles (d= 15 - 40 nm) were grown in samples 
annealed in air and then in molecular hydrogen (samples AH). The largest particles (d= 
40 - 65 nm) were grown in the samples annealed in the H2 (samples H). Based on the 
analysis of TEM images we conclude that the distance between the nanoparticles is vastly 
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larger than the size of the particles, therefore, the electromagnetic interaction between the 
particles is negligible 1 . The particles have a spherical shape and are homogeneous - there is 
no evidence of a core-shell structure of the particles. Annealing in air leads to the oxidation 
of nanoparticles and formation of copper oxide (Cu20)^. However, CU2O nanoparticles have 
not been observed in TEM images, most probably, due to their lower contrast as compared 
to the copper nanoparticles. 

The analysis of TEM results shows that nanoparticles are characterized by Gaussian- 
like size distribution with clear-featured maximum. The dispersion of the distributions is 
quite small. Indeed, for sample Al annealed in air d = (6.1 ± 0.2) nm (e = Ad / < 
d >=0.033), for sample AH1 annealed successively in the air and hydrogen d = (33.9 ± 0.8) 
nm (e = Ad / < d >=0.024), and for sample H 2 annealed in hydrogen d = (51 ± 2) nm 
(e = Ad / < d >=0.039). Because of the narrow size distribution of the nanoparticles its 
influence on the optical response of the nanoparticles can be neglected. 

IV. ABSORPTION OF Cu/Si0 2 NANOCOMPOSITES. RESULTS AND DISCUS- 
SION 

The size-dependence of the width and energy of surface plasmon absorption band at the 
temperatures of 293 K and 77 K we reported earlier in |l9j. In this paper we discuss the 
temperature dependence of SPR absorption bandwidth and its energy. The experimental 
absorption spectra of samples with average diameters of Cu nanoparticles 16.7 nm and 40 
nm) at the different temperatures are shown in figj2j As one can see, the increase in the 
temperature leads to the broadening of the SPR peak. This broadening is well known for 
SPR in metal nanoparticles^, in bulk metals^, and in bulk semimetals (graphite)^.. The 
broadening is due to the increase of the frequency of the electron-phonon scattering (r^) 
with increase of temperature. We have also observed a strong dependence of SPR energy 
on the temperature, see figJ2J 

Fig-El depict the temperature dependences of SPR absorption bandwidth for small parti- 
cles (with diameters smaller than 20 nm), where the gradual (without any features) increase 
of plasmon peak width with temperature takes place. This is similar to the dependences 
observed for SP in small (1.6 - 20 nm) gold nanoparticles^- and SP in bulk silver 22 . The ob- 
served increase of SPR bandwidth for small Cu nanoparticles is close to linear dependence 
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that has been predicted for the temperature dependence of frequency of electron-phonon 
scattering for temperatures higher than the Debye temperature^: Too oc T. However, for 
larger (with the diameters larger than 20 nm) Cu nanoparticles we have observed a gradual 
increase of the bandwidth along with the abrupt significant increase (about 1.5 times) of the 
bandwidth with the increase of the temperature. Such an abrupt increase of the electron- 
phonon scattering frequency was reported in Ref. y at the transition of bulk metals from 
the solid to liquid phase, i.e. at the bulk metal melting. So, we can suppose that such jump- 
like increase of SPR bandwidth is caused by melting (surface melting as we show below) of 
the copper nanoparticles. Similar low-temperature melting of the copper nanorods (at the 



26 



temperatures much lower than copper bulk melting point 1084.5 °C) was observed in Ref. 
by using TEM, SEM and XRD techniques. The ambient environment is an important factor 
affecting the process of nanoparticle melting. In our work we study the copper nanoparticles 
embedded in silica matrix. Therefore, it is reasonably to suppose that the wetting of the 
silica matrix by the molten copper should facilitate the melting that should decrease the 
melting temperature. It is well known that the corresponding abrupt changes must indi- 
cate a first-order phase transition. Let us emphasize that such feature in the temperature 
dependence of SPR bandwidth appears for rather large copper particles and is absent for 
small particles. We note that the temperature of the first-order transition in large copper 
nanoparticles decreases with decrease of the particle size, see figJSl 

The measured temperature dependences of SPR energy (figHJ) are different for the small 
and large Cu nanoparticles similar to the temperature dependence of plasmon bandwidth. 
For small Cu particles the monotonic red shift with the increase of temperature is observed, 



similar to small (1.6-20 nm) gold nanoparticles reported in Ref. [21. For large copper 
nanoparticles the temperature dependence of SP absorption band energy demonstrates the 
non-monotonic character, namely at lower temperatures an absorption band is characterized 
by red shift with the increase of temperature and at higher ones the blue shift is observed. 
Let us note that the temperature (marked as T m ) at which the red shift changes to the blue 
one decreases with decrease of the particle size. 

Let us discuss the temperature dependences of the SPR bandwidth and energy. The 
observed changes in spectral characteristics of SPR (energy and bandwidth) should man- 
ifest some transformation in the structure of nanoparticle. The transformation can be ei- 
ther the solid-liquid phase transition (melting) or the solid-solid phase transition (struc- 
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tural transition). The later was observed by electron diffraction technique for free copper 
nanoparticles^ 1 ^ and by HRTEM technique for free gold nanoparticles^. First, let us discuss 
the melting as possible reason for the observed temperature dependences of SPR spectral 
characteristics. As we mentioned earlier, the jump-like increase of the SPR bandwidth at 
certain temperatures suggests the first-order phase transition solid-liquid (melting) in copper 
nanoparticles with diameters larger than 20 nm. Above dependence is quite similar to the 
dependences observed at the melting of bulk metals^^i. The melting of such kind begins 
by the process of the melting of thin layer close to the surface, i.e. surface melting. The 
surface melting generally occurs at temperatures much lower than the bulk melting point. 
A stable coexistence of spatially separated solid and liquid phase takes place during the 
increase of temperature up to the melting point. Melting of entire volume of the particle 
occurs at the volume melting temperature. The surface melting has to manifest itself in the 
SPR bandwidth mainly. The volume melting has to affect mainly on the energy of SPR as in 
copper the dipolar mode of SP resonance coincides with the dominant interband transition 
edge. One of the important characteristics of the bulk melted metals is the blue shift of the 
SP band occurring at the increase of the temperature^. That effect is due to the reduction 
of polarizability of the melted ion core. In solid state the red shift of SPR energy is ob- 
served due to lattice dilation occurring at the increase of temperature. The lattice dilation 
decreases the electron density, and, hence, the plasma frequency uj p = n e e 2 /eom e . Therefore, 
one can conclude that at the temperatures higher than temperature at which the red shift of 
SPR energy with the temperature changes to blue one the copper nanoparticles are melted, 
i.e. this temperature (marked as T m ) is the melting temperature of the Cu nanoparticle. 
Correspondingly, the temperature of jump-like increase of SP bandwidth is the temperature 
of surface melting (marked as T sm ). At the temperatures lower than T srn the nanoparticle 
is solid entirely. In the range T sm <T< T m the coexistence of solid core and liquid surface 
occurs. The dependence of the temperatures of melting and surface melting on the diameter 
of Cu nanoparticle is shown in fig. [5j Let us note that these temperatures are equal for 
smaller particles with diameters up to 40 nm, and for larger particles the melting tempera- 
ture exceeds the surface melting one. Most probably, that is due to greater role of surface 
in small particles. 

Let us discuss the solid-solid phase transition as the possible cause of the obse rved t em- 
perature dependences of SPR spectral characteristics. As it was shown in Ref. 



27 



2M the 



s 



morphology of small copper nanoparticles (smaller than about 4 nm) is icosahedral (Ih) or 
decahedral (Dh), where Ih and Dh are noncrystalline structures. The morphology of the 
particles larger than 4 nm can be Ih, Dh or fee (face centered cubic), where fee structure is 
crystalline structure, i.e. the structure of bulk copper. The similar situation occurs in gold 
nanoparticles 29 where the morphologies of small nanoparticles (smaller than about 6 nm) 
are noncrystalline, namely icosahedral (Ih) or decahedral (Dh). The morphology of gold 
particles larger than 6nm can be Ih, Dh or fee (face centered cubic), where fee structure 
is crystalline structure of gold. The Ih and Dh structures are noncrystalline, i.e. these 
structures do not exist for bulk crystals. Since the copper nanoparticles can exist in form 
of three different structures (Ih, Dh and fee), there is a possibility, that the observed pecu- 
liarities in temperature dependence of SPR bandwidth and energy ctt T sm cincl T m are the 
results of some solid-solid structure transition, i.e. Ih-Dh, Ih-fcc or Dh-fcc. However, the 
solid-solid transitions between noncrystalline Ih and Dh phases as well as transition from 
noncrystalline Ih or Dh phases to crystalline fee one seems to be of low possibility due to the 



following reason; as it was shown in Refs.27.28.29 the relative fractions of the particles with 



noncrystalline Ih or Dh structures are vanishing for large particles with diameters higher 
than 20 nm, where the above mentioned peculiarities in temperature dependences of SPR 
energy and bandwidth are observed. For sizes higher than 20 nm nearly all the particles 
have the crystalline fee structure. Therefore, we believe that the melting is a more probable 
reason for jump-like increase of SPR bandwidth at T sm and nonmonotonic behaviour of 
SPR energy at T m observed for large Cu nanoparticles than any kind of solid-solid struc- 
tural transition. Since the temperature behavior of the energy and bandwidth of SPR at the 
melting of copper nanoparticles larger than 20 nm is the same as for bulk copper, one can 



es_have crystalline fee morphology, 
291 . Their melting occurs bulk-like, 



assume that, similar to bulk copper, the large Cu partic 
which is in full agreement with the results of Refsj27 
i.e. it is first-order phase transition. 

Thus, an abrupt jump-like increase of SPR bandwidth at T sm has to be attributed to the 
surface melting of studied large Cu nanoparticles. It is noted above that the temperature 
dependences of energy and bandwidth of SPR for small copper nanoparticles are quite dif- 
ferent from ones for large particles. Namely, for small particles an increase of the bandwidth 
with temperature is gradual and is not characterized by abrupt jump at certain temper- 
ature, and the temperature dependence of SPR energy is monotonically decreasing. Both 
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these dependences are similar to the corresponding dependences observed e.g. for small gold 
nanoparticles^-. Let us note that temperatures reached in experiments with the small Cu 
nanoparticles are the same as temperatures at which evidences for large Cu nanoparticles 
melting have been obtained. So, as the melting temperature is lower for smaller particles, one 
can suppose that the melting of small nanoparticles occurs in other way than for large ones. 
The molecular dynamics simulations of the solid-liquid phase transitions in various metal 
nanoparticles were performed in Refs. where it was shown that two characteristic 

particle sizes exist. For particles with size larger than first one (t^) the solid-liquid phase 
coexistence is stable. The melting of such particles is first-order transition, i.e. the melting 
has to occur likely to bulk metals. For particles with size smaller than second characteristic 
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that such small 



value (d\) the solid-liquid coexistence is unstable. It is shown in Ref. 
particles can be superheated and non-melted up to temperatures exceeding the bulk melting 
temperature. Finally, in the middle-sized particles the coexistence of the solid and liquid 
phases is met ast able, and dynamic and transient coexistence can exist. Taking into account 



the results of Refs. 



27 



28 



29l . one can suppose the origin of the critical size of about 20 nm 



obtained for our copper nanoparticles. As it was shown in these studies, the relative fraction 
of both copper and gold nanoparticles with fee crystalline morphology increases with the 
increase of the particle size, and on the contrary the relative fractions of the noncrystalline 
Ih and Dh morphologies decrease. In the region of size higher than 20 nm the nanoparticles 
with Ih and Dh structures are absent, and the nanoparticles with fee structure only exist. 
Therefore, it is plausible to assume that, since the copper nanoparticles larger than about 
20 nm have the crystalline fee morphology, their melting would occur bulk-like. The copper 
particles smaller than about 20 nm have mainly the noncrystalline Ih or Dh morphologies. 
Therefore, their melting occurs in other way than in bulk-like large particles. So, one can 
suppose that the size of about 20 nm is the c?2 characteristic size mentioned above. 



Therefore, taking into account our experimental data, data of Ref. 
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28 



J29l and data of 

computer simulations^ 1 ^ one can conclude that the melting of large studied Cu nanoparti- 
cles (with diameters larger than 20 nm) is first-order transition, and correspondingly, it is 
similar to the melting of bulk copper, i.e. the stable solid-liquid phase coexistence occurs 
in the temperature range from surface melting to volume melting. And, most probably, Cu 
nanoparticles with diameters up to 20 nm are characterized by the unstable or metastable 
solid-liquid phase coexistence. So, the observed monotonic temperature dependences of the 
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surface plasmon bandwidth and energy suggest either the existence of the solid superheated 
particles or the fact that melting of such particles occurs gradually, and it is not first-order 
transition, that is quite different from the melting of bulk copper. We note, that the results of 
molecular dynamics that we use in our analysis of the observed phase changes are conducted 
for isolated nanoparticles whereas our Cu particles are embedded in a matrix. In principle, 
the behavior of isolated nanoparticles can differ to embedded particles. However, we think 
that environment (silica matrix) can change the values of melting characteristics (critical 
size of nanoparticle, melting temperature), but the qualitative behavior of the nanoparticle 
at melting should remain the same. 

V. CONCLUSIONS 

In conclusion, we studied the temperature dependences of the bandwidth and energy of 
surface plasmon in spherical copper nanoparticles embedded in silica matrix. The obtained 
dependences indicate that for small particles (smaller than 20 nm) either the solid super- 
heated nanoparticles exist or the melting is not a first-order phase transition that is different 
from the melting of bulk copper. The melting of large Cu particles (larger than 20 nm) is a 
first-order transition accompanying by jump-like increase of surface plasmon bandwidth and 
non-monotonic dependence of surface plasmon energy with the temperature that is similar 
to bulk copper melting. 
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FIG. 1: TEM images of Cu nanoparticles in silica matrix, a - Al sample annealed in air (mean 
diameter of copper nanoparticles < d > = 6.1 nm), b - AH2 sample annealed successively in air 
and hydrogen (< d > = 16.7 nm), c - HI sample annealed in hydrogen (< d > = 50.5 nm). 

29 K. Koga, T. Ikeshoji, and K. Sugawara, Phys. Rev. Lett. 92, 115507 (2004). 
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FIG. 2: Absorption spectra of the copper nanoparticles with mean diameter of 16.7 nm (a) and 
nm (b) in silica matrix measured at different temperatures. 
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FIG. 3: Temperature dependence of the surface plasmon absorption bandwidth (normalized to SP 
bandwidth at 77 K) for Cu nanoparticles of different sizes. Arrows mark the temperatures of the 
surface melting of nanoparticles. 
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FIG. 4: Temperature shift of the SP absorption band from the spectral position corresponding to 
77 K for copper nanoparticles of different sizes. Arrows mark the temperatures of the melting of 
nanoparticles. 
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FIG. 5: Dependences of the temperatures of melting and surface melting of copper nanoparticles 
on their diameter. 
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